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We develop a microscopic model for the recently demonstrated double quantum dot (DQD) maser. 

In characterizing the gain of this device we find that, in addition to the direct stimulated emission 
of photons, there is a large contribution from the simultaneous emission of a photon and a phonon, 
i.e., the phonon sideband. We show that this phonon-assisted gain typically dominates the overall 
gain which leads to masing. Recent experimental data are well fit with our model. 
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The coherent generation of light in a laser provides 
fundamental insights into the interaction between light 
and matter [Tj . Lasers operating in the few-emitter limit 
probe this interaction at the level where quantum effects 
are crucial for understanding the device operation [mi- 
single emitter lasers were first demonstrated in atomic 
systems OE] and, subsequently, extended to solid-state 
systems, where one must contend with a strong coupling 
of the emitter to the surrounding environment [ZHIQ]. 

Several groups have explored the possibility of achiev¬ 
ing a maser with gate-defined semiconductor quantum 
dots as the gain medium diHiii- Recently a successful 
demonstration of such a maser was achieved by coupling 
two InAs nanowire double quantum dots (DQDs) to a 
microwave cavity [15]. Due to the large Coulomb charg¬ 
ing energy Ec ^ b meV, these systems provide tunable 
gain from GHz to THz frequencies using external gate 
voltages. Operating in the few-emitter limit, they may 
enable the creation of quantum states of light m and 
entangled states of DQDs and light [T7H2Q] . The strong 
environmental coupling in these devices allows the study 
of competing emission mechanisms, e.g., phonon versus 
photon [21]. While the role of electron-phonon coupling 
has been considered in previous work on optical quantum 
dot lasers [221125] , electrically driven quantum dots probe 
a much lower energy scale. Finally, previous theoretical 
work predicts a small, narrow gain feature in the DQD 
emission spectrum mHE]. This is in contrast with the 
experimental results, where high gain is observed over a 
much larger energy range [HITS]. Resolving this discrep¬ 
ancy is crucial for future applications of the DQD-cavity 
system to both maser operation and quantum informa¬ 
tion tasks. 

In this Letter, we develop a microscopic model for the 
recently demonstrated DQD maser m- In characteriz¬ 
ing the gain of this device, we find, in addition to the 
direct stimulated emission of photons into the cavity, a 
large contribution from transitions that involve the si¬ 
multaneous emission of a photon and a phonon, i.e., the 


phonon sideband. These effects have not been consid¬ 
ered in previous related work [TI11I3] . Under typical ex¬ 
perimental conditions, the phonon sideband dominates 
the gain and, therefore, sets the energy range over which 
masing occurs. We find the experimental data from Ref. 
na are well fit with a theoretical model accounting for 
this phonon process. 

A schematic of a DQD maser is shown in Fig. 1. The 
gain medium consists of one or several DQDs coupled to 
the common mode of a microwave resonator [Fig. 1(a)]. 
With a bias applied across the DQDs, current flows via 
single electron tunneling and, in Ref. mmsi , gain was 
observed in the cavity transmission. However, this gain 
occurred over a much wider range of DQD transition fre¬ 
quencies than the cavity resonance and was much larger 
than predicted from a Jaynes-Gummings model. We can 
understand the broadening of the gain at a qualitative 
level by noting that the electron-phonon interaction will 
dress the electronic states of the quantum dot with the 
phonons in the nanowire. This leads to a phonon side¬ 
band whereby energy is conserved through the simulta¬ 
neous emission of a phonon and cavity photon, as illus¬ 
trated in Fig. l(b-c). To understand how this affects 
the maser note that effective maser operation requires a 
large photon emission rate, a large population inversion, 
and a rapid repumping rate. The peak emission rate for 
the direct process occurs when the DQD is on resonance 
with the cavity. Without precise tuning of the system, 
this does not always correspond to the optimal operating 
point for the maser (e.g., as is the case in Ref. [HdS]). 
Furthermore, in the presence of charge noise it is diffi¬ 
cult to stabilize the DQD at the resonance condition. In 
this far-off resonant regime, we show that the phonon 
sideband strongly dominates over the gain from direct 
photon emission. As a result, the DQD maser dynamics 
is typically dominated by this phonon-assisted process 
and not direct photon emission. 

DQD Gain Medium - To analyze the masing process 
we first need to characterize the DQD gain medium and 
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FIG. 1: (a) Schematic of combined DQD, cavity, and phonon 
system, (b) DQD energy spectrum versus detuning e for tc = 
16.4 /ieV, the direct photon emission process, and the phonon 
sideband for — 8 GHz. (c) Transport cycle for DQD 

maser: a hnite source drain bias leads to current flow via 
single electron tunneling. The inter dot charge transition is 
accompanied by direct photon emission into the cavity (near 
zero detuning) or a second order process involving emission 
of a photon and a phonon (large detuning). Electrons tunnel 
onto(off of) the left(right) dots at rate r^(r)- 


its coupling to the microwave cavity. Following previous 
theoretical work HIHIS], we develop a simplified micro¬ 
scopic model for the system, consisting of a DQD, single 
mode cavity, phonon bath, and leads. From this model 
we can extract the gain, which is determined by three dis¬ 
tinct quantities: the DQD photon emission and absorp¬ 
tion rates and the population inversion of the DQD. We 
then extend this result to multiple dots in the sideband- 
dominated limit. 

Due to Coulomb blockade, each DQD can be restricted 
to two orbital states \L) and |i?), where \L) has (M + 
l^N) electrons and |i?) has {M^N + 1) electrons in the 
(left, right) dots. The different charge configurations of 
these states results in an electric dipole moment on the 
order of D ^ 1000 eao, where e is the electronic charge 
and ao is the Bohr radius. The Hamiltonian describing a 
single DQD coupled to a cavity is given by 

Ho = ^cTz + tc<7x + hixjcO'a + hgcaz{a + a^), (1) 

where are Pauli-matrices operating in the orbital 
subspace \L) and |i?),e is the detuning between the 
two dots, tc/h is the interdot tunneling rate, ujc is the 
cavity frequency, Qc is the DQD-cavity coupling, and 
ai(a) are cavity photon creation(annihilation) operators, 

. The electron-phonon interaction takes the generic form 
Hep/h = + dqu), where 

uJu{q) is the phonon dispersion, Xu{q) is a coupling con¬ 
stant that depends on momentum q and mode index u 
and are phonon creation(annihilation) opera¬ 

tors. The exact form of \y{q) is set by the electronic 
wavefunctions, material properties, and boundary condi¬ 
tions. We focus on the phonon properties of nanowire 
QDs [26] . 

Diagonalizing the first two terms in Hq leads to the 


eigenstates |±) 

1+) = cos(6>/2)|L) - sin(6>/2)|i?), (2) 

|-) = sin(6>/2)|L) + cos(6>/2)|i?), (3) 

where 0 = tan“^(2tc/e). These states have an energy 
splitting fujJd = shown in Fig. 1(b). Writing 

the Pauli-matrices in this new basis, the interaction be¬ 
tween the DQD, phonons, and cavity photons is 

Hint = h{cos 9az + sin6'(73;) [fife a + El ^v{q)aqv + h.c]. 

q,u 

From this interaction we see that both phonons and pho¬ 
tons will cause relaxation from |+) to |—); therefore, sin¬ 
gle electron tunneling through the dots will be correlated 
with photon and phonon emission [26l |27] . 

In the presence of a finite source-drain bias, an electron 
first tunnels from the drain to the left dot, followed by 
an interdot charge transition from \L) to |i?), and then 
leaves the right dot by tunneling to the source. In the 
context of the maser, this can lead to population inver¬ 
sion when e > 0 as it continually repumps |+). In the 
limit where only single electrons can tunnel through the 
DQD, this process can be modeled by including a third, 
empty dot state |0) with incoherent tunneling rates F^ 
from |0) ^ \L) and F^ from \R) |0) [see Fig. 1(c)]. 

Thus the dynamics for a single DQD can be described by 
the master equation for the density matrix p [12] 

f> = dH,p] + KV[a]p + r, V[\L)mp + r, V[\0){R\]p, 

where H = Hq Hep describes the coherent dynam¬ 
ics (including the phonons) and the incoherent evolution 
is described by the Lindblad super-operators VlAjp = 
— 1/2{A^A, p} + ApA\ for any operator A, correspond¬ 
ing to cavity decay, at rate k, and inelastic electron tun¬ 
neling. 

Neglecting the phonons, the emission rate of photons 
into the cavity can be found perturbatively for small pc by 
using the Heisenberg-Langevin equations for the DQD- 
cavity system in a rotating wave approximation: 

(7_ = -(F + i(u;d -cJc))cr- -h ipc sinO a ^cFz^d, (4) 
a =—t<i/2a ^ igcAnO (7-^ He- (5) 

Here J-c{d) the associated noise operators for the cav- 
ity(dot) baths and F is the total dephasing rate (de¬ 
fined below). Adiabatic elimination and mean field 
theory, i.e., (acr^) ~ (a) {(Tz)^ appropriate for large F, 
gives the equation of motion for the cavity photon num¬ 
ber Tie =< a^a >, Tie = —(/^ — R {(7z)){a^a). Where the 
direct photon emission rate for the DQD is 

gc p 

^ 0J^^T^ + {cod - • 


( 6 ) 
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The dominant effect of phonons is to induce relaxation 
from 1+) to |—) via phonon emission. Neglecting cavity 
effects, the zero-temperature emission rate is given by 
Fermi’s Golden rule as 

7. = (7) 

where — cj) is the spec¬ 

tral density of the phonons. In the presence of ther¬ 
mal phonons with distribution np{uj), the total emission 
rate is 7 ^ = + !)• Starting in |—), there is 

also absorption at the rate 7 ^ = 7 d%(<^d)- Using these 
expressions, we can write the total dephasing rate as 
^ ~ + 7 I + rr)/ 2 . We treat the phonon spectral 

density J{uj) using the microscopic model and measure¬ 
ments in Refs. [26l|28]. We take a 25 nm radius nanowire 
with a separation between the two dots d = 120 nm, an 
axial confinement a = 25 nm for each dot, and a phonon 
speed of sound = 4000 m/s [29] . 

The effect of phonons on the photon emission is cal¬ 
culated by first performing a polaron transformation 
H' = UHW with [28] 

JJ _ ^[9cia-a^)/uJc+J2q,^ cosOa^ ^ 


which removes the az terms in the interaction. The po¬ 
laron transformation serves to dress the electronic states 
of the DQD with the ambient phonons in the environ¬ 
ment. Perturbatively in gc/(^c and \y{q) /ujy{q) this re¬ 
sults in explicit terms in the Hamiltonian, which have 
not been considered in previous theoretical treatments of 
the DQD maser [TTHT3] . describing second order photon- 
phonon processes 


, _ 4tc6 igcXi^iq) [//>,, w t t \ 

+ (w^(g) - u:c){a a\^ - a'aq^) ay + H", (9) 


where H" contains terms that do not directly couple 
photons and phonons. The first term in Eq. ^ leads 
to phonon assisted emission [Fig. 1(b)], whereby relax¬ 
ation from 1 +) to |—) occurs by emitting a phonon of 
frequency ujy{q) = along with a cavity photon. 

The second term leads to phonon assisted absorption, 
whereby relaxation occurs by emitting a phonon of fre¬ 
quency ujjy{q) = ujd ^ ujc and absorbing a cavity pho¬ 
ton. Using Fermi’s golden rule, these two terms give 
the zero-temperature, phonon-assisted-photon-emission 
7 e and absorption 7 ^ rates as 


7e 


7a 


Ojj ujI {oJd - UJc, 


2 ^ ^0)5 




J {^d + ^c)- 


( 10 ) 

( 11 ) 



FIG. 2: (a) Population inversion (cr^), direct photon rate 

R, and phonon-assisted emission 7e and absorption 7a rates 
as a function of e with tc = 16.4 /icV, gd^'^ — 90 MHz, 
Ct;c/27r = 8 GHz, J(2tdK) — 5 GHz, = 4 GHz, 

and Teff = 0. (b) 7e/^ plotted as a function e and The 

dashed line corresponds to — 8 GHz used in Ref. [15] . 


In the presence of thermal phonons, we define 7 ^’^ anal¬ 
ogously to case for the direct phonon process, where ^ 
refers to transitions from |+) to |—) and vice-versa for 
t- These thermal contributions are important because, 
in addition to the ambient thermal phonons in the DQD, 
pumping current through the dot will generate a large 
population of phonons through Ohmic heating of the 
nanowire. Since the equilibration time of the phonons 
is on the order of picoseconds a/c^ ^ 10 ps) and the 
cavity dynamics occur over a timescale of hundreds of 
nanoseconds ^ 100 ns), we can take the phonon 

bath to be in equilibrium with an effective temperature 
Teff, such that np{uj) = _ i)-i. In Kef. [T4] . 

Teff was estimated to be as high as 1 K due to the large 
nA currents flowing through the nanowire. 

Figure 2(a) shows the key quantities in determining the 
gain for parameters similar to Ref. [HUS]: the popula¬ 
tion inversion {(Tz)^ obtained from Eq. (5) in the absence 
of the cavity, and the various photon emission and ab¬ 
sorption rates. For e = 0 the DQD eigenstates are equal 
admixtures of \L) and |T) and the inversion is small, 
while for large, positive e, |-h) ^ |T), the system be¬ 
comes completely inverted as seen in Fig. 2(a). Although 
R and 7 e are comparable in magnitude, R dominates at 
small detunings, where the population inversion is small, 
and 7 e dominates at large detunings, where there is large 
population inversion. Based on Eq. ( 6 |[lQ), it is possible 
for R to dominate the gain at large e when ujc ^ 2tc. 
However, this analysis has so far neglected charge noise 
in the system. In Ref. m IIS, this was estimated to 
lead to slowly varying noise in e with an rms value of 
(20 — 40) /ieV ^ hT ^ 1 geV. In the presence of such 
large noise, 7 e will dominate over R for large e. Fig¬ 
ure 2(b) shows 7 e for varying e and cJc- Because the 
phonon-assisted process is perturbative in gd^a H has 
the strongest effect for small cavity frequencies. The sec¬ 
ond peak at e = 150 yueV arises from the second phonon 
branch in the nanowire. 

DQD Maser - Away from the masing threshold, we can 
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FIG. 3: (a) (Circles) Experimentally measured gain in 

one of the DQDs with — 8 GHz, /^/27r = 2.6 MHz, 

4 = 50(10) //eV, r^/27r = r^/27r = 17(2) GHz, Teff = 3(1) K, 
and ^c/27r = 100(20) MHz. (Black) Fit to the theory with 
J(2tc/^) = 2.4(2) GHz a free parameter and only including 
contributions from the lowest phonon branch. (Blue) Includes 
contributions to J{uj) from the second phonon mode and sub¬ 
strate phonons m- (Dashed) Gain neglecting the phonon- 
assisted contributions, (b) Gain rate function g{uo) plotted 
versus e and with parameters as in Fig. 2(a). The sec¬ 

ond peak near e = 150 /xeV arises from the second phonon 
mode in the nanowire. 


find the response of the system within mean field theory. 
Including thermal effects, the Heisenberg-Langevin equa¬ 
tions gives rise to the mean field equations for the field 
amplitude a = (a) and the population in the upper state 

= (l+)(+l) HSIEI]: 

a =— g{u)\/2 ^ id)a ^ ( 12 ) 

u = rp(i4o — u) — S{u) \a\^ , (13) 

where we have defined the gain rate function g{u) and a 
saturation function S{u) as 

g = R{2u-l) + {Y - 71)m - {Ye - 7l)(l - u), (14) 
S = R{2u-l) + {Y + Y)u - {Ye + 7l)(l - «)• (15) 


Here we have introduced the drive with amplitude Vi and 
frequency uji. The detuning 5 = Uc — ooi — RV(2u — 
1)/2A, includes the cavity line pulling [T]. F^ and uq are 
the effective pumping rate and upper state population, 
respectively. The full expressions are given in Ref. m 
[13]. For large e and F^^^, they reduce to Tp « (F^ + 

2 F^) and ~ 1 “ 2 ( 7 ^ + ld)/^P' 
experiment, where there are multiple DQD (two), the 
large dephasing rate F allows u to simply be replaced by 
the average upper state population in each DQD and g 
to be multiplied by the number of DQDs. 

For weak driving fields and below threshold operation, 
the normalized gain |(a(^; g)f / |<a(0; 0)|^ is given by 


G{6) = 


[k - g{uo)Y + 


(16) 


From transmission measurements it is known u;c/27r = 
8 GHz and nl2ii = 2.6 MHz [T5|, modeling the current 


through the dot at finite bias gives tc = 50(10) /icV, 
F^/27r = F^/27r = 17(2) GHz, Teff = 3(1) K [T4], and 
the gain at zero bias gives 5fc/27r = 100(20) MHz [Ml 
[30] . To account for charge noise, we convolved the gain 
with a Gaussian of width 40(10) jieV [18]. Finally, we 
find J{2tc/h) = 2.4(2) GHz by fitting the gain at finite 
bias including only the first order phonon branch in J{uj) 
(which is a valid approximation for |e| < 200 fieV for a 25 
nm radius nanowire [26]). To match the broad tails in the 
gain data for e > 200 /ieV [see Fig. 3(a)], we include the 
coupling to the second longitudinal mode and substrate 
phonons [29] . 

Figure 3(a) shows the comparison between the mea¬ 
sured G(0) for a single DQD and a fit to our model. From 
the data we can conclusively rule out a model with just 
the direct photon emission process as it would require a 
DQD-cavity coupling gc 10-100 times larger than what 
was measured. On the other hand, when the phonon- 
assisted processes are included, we find good agreement. 

Equations (12)-(13) predict a masing transition when 
g{uo) > n. This is consistent with the experimental re¬ 
sults, where there are two DQDs in the cavity, each with 
peak gain rates slightly below the cavity linewidth. When 
only one DQD is configured to maximum gain, no signif¬ 
icant photon emission is observed; however, when both 
are tuned to maximum gain, such that the combined gain 
rate is greater than the cavity linewidth, masing is ob¬ 
served [T5] . 

Because of the strong dependence of the gain on the 
phonon-assisted process, measuring the gain near thresh¬ 
old is a sensitive measurement of the phonon spectral 
density In particular, by tuning ujc and measuring 

the gain curves as in Fig. 3(a), one could precisely de¬ 
termine the frequency dependence of J{uj) by extracting 
g{uQ). This is illustrated in Fig. 3(b), which shows g{u{)) 
at Teff = 0 for varying e and uOc^ where we see that there 
will be a second peak in the gain at low frequencies when 
equals the gap to the second longitudinal phonon 
mode of the nanowire [29] . 

Finally, this work shows that phonons will be impor¬ 
tant for circuit quantum electrodynamics (QED) experi¬ 
ments involving spin-photon entanglement and the gen¬ 
eration of non-classical states of light. In conventional 
cavity QED, the fidelity of these operations is limited 
by the largeness of the Purcell factor g‘1/ 1 ^ r m- In the 
case of the DQD, as gc approaches uJc the phonon-assisted 
processes can dominate over the bare relaxation rate F. 
This will ultimately constrain the fidelity of these op¬ 
erations, but it also represents an unexplored regime of 
cavity QED that is unique to the solid-state environment 
and energy scales of the DQD system. 
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Supplemental Material 


Heisenberg-Langevin Equations - We give an elemen¬ 
tary derivation of the Heisenberg-Langevin equations 
used in the manuscript, which follows standard treat¬ 
ments found in, e.g.. Refs. uniET]. Neglecting the leads, 
we break the total Hamiltonian following the polaron 
transformation into three separate terms describing the 
DQD-cavity system, the phonon and photon reservoirs, 
and the coupling between the system and the reservoirs 
(fi = 1) 

H = Hs + Hr + Hsr, (SI) 

Hs = + gcSm9{a^a- + h.c.), (S2) 

Hr = (S3) 

q,u n 

Hsr = —R\ a — R\(T- — R\ aV_ — R\ acr_ — h.c. 

(S4) 

where we introduced a bath of modes that couple to the 
cavity field hn and neglected counter-rotating terms and 
higher order terms in (^fc, \,v) in Hs and Hr. Excluding 
Rc^ we implicitly derived expressions for the reservoir 
operators in the main text. Writing them out explicitly 


Rc — ^ ^ : 


(jOfj ^' 

^ q,u 


'qy-) 


_ 4tc6 9c^y{(l) / / N , N 

~ \ P ^c)ciqy, 

q^y ^c^y[Q) 

_ Atce gcKiq) / / N X 

~ ^c)Ciqy‘ 

d ^ ^c^y[Q) 


(55) 

(56) 

(57) 

(58) 


In a frame rotating at frequency u, we can formally ex¬ 
press the reservoir operators in terms of the system op¬ 
erators: 


Rc{t) = 

n 

+ i Tdi' (S9) 

'^0 n 
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and similarly for the other R^. The key approximation 
is to assume that the system operators change slowly in 
time so that they can be removed from the integral in 
the second term of Eq. ( |S9| ) . This allows us to write 


Rc{t) « 

n 

K = 2TT'^\Tnf S{0J - W„). 


(510) 

(511) 

(512) 


and similarly for the other R^. With these expressions, 
the equations of motion for the system operators close 
and we can solve for the dynamics solely in terms of the 
system operators (provided we know the correlations of 
the For example, the Heisenberg-Langevin equation 
for cr_ takes the form 


&- = —(T + iA)cr_ + igc sinO aaz + cTzRd (S13) 

- cr^a(YaV_ - - cr2a'^(ya(7_ - Ja), 

where A = ood — oOc is the detuning of the dot from the 
cavity (including corrections from the polaron transfor¬ 
mation) and Re{a) are the noise operators associated with 
the phonon-assisted emission(absorption) terms. They 
satisfy (J’t (t) Defining 7 ;^(f) = 

f* we can express CT-(t) to second 

order in gc as 

ioc sin 0 aaz , . + . . , . 

O’- ~ „ , .. -h aa^rieit) + a'a^^rjait) + a^rid{t). 

L IA 

(S14) 

Inserting this into the Heisenberg-Langevin equations for 
a and cr^, gives the mean field equations for the field 
amplitude a = (a) and the population in the upper state 
u = (I+)(+!) shown in the main text. 

Phonon Spectral Density - Here we discuss the fre¬ 
quency dependence of the phonon spectral density, in¬ 
cluding higher order modes and coupling to substrate 
phonons. For the nanowire phonons, we can express the 
matrix elements for the electron-phonon interaction as 
[261132] 


= {L\V^AL) - {R\V^,\R) (S15) 

= Jdx{\(f>e{x)f - \(f>r{x)f)Vqi,{x)e^‘>^ 

= 2isin{qdf2)e~'^^^^‘^Mqjj, 

Mg, = jdx \c^{x)\^ Vg,{x)d'>\ (S16) 

where Vqj, is the interaction potential for the phonon 
mode with quasimomentum q along branch z is the di¬ 
rection along the nanowire axis, d is the spacing between 
the dots, and (pi^r) = ^ Ihe envelope wavefunctions 
for the electrons in the dots (assumed to be identical). 


We approximate the envelopes by a cylindrically sym¬ 
metric Gaussian 0(cc) = (ay^)”^/^ (/)^(r) with 

a = 25 nm. The interaction potential can be broken up 
into a contribution from the deformation potential 
and the piezoelectric potential where the deforma¬ 
tion potential is given by [33] 

= -D_lV± • Ug,{x) + iqDz z ■ Uq,{x), (S17) 

where Di are the deformation potential constants and 
Uqjy is the displacement of the phonon mode such that 
the quantized displacement field is given by Uu = 
^- 1/2 Uqj^e^^^ [aqy + 0^-qy) {N is the number of unit 
cells). The piezoelectric potential is more difficult to cal¬ 
culate; however, for small q it only depends on the trans¬ 
verse variables and is independent of q. Expanding Vqy 
and performing the integral over space we find 


“ Y 2Muj,{q) 

where M is the average mass of the unit cell and /3y and 
Ey are deformation and piezoelectic constants, respec¬ 
tively, which depend on the different phonon branches u 
and the envelope wave function. In the frequency range 
relevant to the experiment, a rough estimate shows that 
the piezoelectric term should dominate over the deforma¬ 
tion potential term [26] . 

The phonon spectral density for the nanowire is given 
by 

J{uj) = 4'^sm^{qd/2) \Mq,f S{uj - Wq,) (S19) 

q,u 


For a 25 nm radius nanowire, the higher energy waveg¬ 
uide modes of the nanowire have a gap of ^ 200 jaeV/h 
so there is only a single phonon branch for energies be¬ 
low this cutoff [26]. Furthermore, for small q this mode 
is longitudinally polarized and has a linear dispersion 
which implies 


J(cu) ~ Jj2 


sm^(u)dl2cn) ^_^2^2/2^2 

Cudj2Cn 


(S20) 


In the data in Fig. 3(a) of the main text, there are 
broad tails in the gain for large detuning. At such large 
frequencies our approximation for J(cj) including only 
the first order mode needs to be extended to include 
higher order modes. We include the second order mode 
by taking a dispersion of the form 002 (d) = 
which gives rise to the contribution 

(S21) 

q 2 \pj)(i/ A 

where 0 is the Heaviside step function and g 2 (<^) = 
up — uJq. In addition, the phonons in the amor¬ 
phous silicon nitride (SiN^,) substrate will begin to hy¬ 
bridize with the nanowire phonons, as was noted in Ref. 
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[26] . This leads to a contribution from the piezoelectric 
potential [28] 

Jsp{^) ^ (l — sinc(co’(i/cs)), (S22) 

which we expect to dominant over the deformation po¬ 
tential for the considered frequencies. Here Cs ~ 11000 
m/s is the longitudinal phonon speed of sound in SiN 


and sinc(x) = sm{x)/x. Both these terms give a large 
contribution to the spectral density, whereas the contri¬ 
bution from the first order mode becomes exponentially 
suppressed at these frequencies. Adding these terms to 
J{uj) with hwo = 200 yueV, we are able to account for the 
broad tails in the data in Fig. 3(a). 


